Pairing Symmetry and Spin Excitations of Iron Selenide Superconductors 
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We propose an effective single-band model for the newly discovered iron selenide superconductors 
AjFe2_vSe2 (A=Tl,K,Rb,Cs). Based on this minimum model and the random phase approximation, the d,2_,:2 
pairing symmetry is revealed theoretically, which may be understood in the framework of Fermi surface topol- 
ogy. A common origin of superconductivity is elucidated for this compound and other high-Tc materials. The 
spin excitations at (n, n/2) in superconducting states are observed, in good agreement with the neutron scatter- 
ing experiments. The spin resonance is indicated to show up only for the d-wave pairing, which provides an 
additional indication for the rf-wave pairing symmetry in this family of superconductors. 

PACS numbers: 74.70.-b, 74.20.Rp, 74.25.Ha 



I. INTRODUCTION 



Recently, the discovery of a new type of iron-based super- 
conducting (SC) compounds AiFe2-vSe2 (A=Tl,K,Rb,Cs) has 
attracted much attention 1 1-3]. Both the angle-resolved pho- 
toemission spectroscopy (ARPES) experiments ||4Mland the 
local density approximation (LDA) calculations ll9l- [llll have 
reported that the hole Fermi surface (FS) pockets disappear 
and only electron ones exist. Thus this compound may be 
heavily electron-doped, which challenges the current physi- 
cal picture for iron-based SC materials. Due to the absence 
of the hole-like FS, obviously, the inter-band FS nesting pic- 
ture 1 12*, Ts'], which seems to work well for iron pnictides, is 
unable to straightforwardly account for the superconductivity 
of this compound. The pairing symmetry is also an impor- 
tant and critical issue to be addressed, which may be different 
from that of iron pnictides due to the absence of the hole FS 
pockets. Actually, the ARPES experiments ||4)-(7t] indicated 
that the FS pockets around the M = (;r, 7r)-point are fully 
gapped while near the F - (0, 0) point the gap is negligibly 
small |4, 5], which is contrast to the i,2,,:-wave symmetry. Up 
to now, the pairing symmetry is still unclear: both i-wave and 
cZ-wave symmetries have been proposed 

Due to the presence of magnetism in the parent compound, 
usually the spin excitation is proposed to account for super- 
conductivity in high-Tc materials lfl2l [Tsll . It is natural and 
important to ask whether this picture still works well in iron 
selenides. This issue was first studied theoretically based on 
the five-orbital model. The d-wa\e pairing symmetry, which 
meditated by the spin dynamics, was proposed 1 14i, J5il : while 
the results of spin susceptibility are different, with the (tt, n)- 
spin excitation proposed in Ref. 111411 and {n, 0.6257r)-one pro- 
posed in Ref. [ 15.1 ■ 

Experimentally the inelastic neutron scattering (INS) ex- 
periment is one of the most powerful ways to study the spin 
excitations and can provide us directly the imaginary part of 
the spin susceptibility in the whole momentum and frequency 
space. The INS experiments on iron selenides have been re- 
ported by several groups. For the insulating parent compound, 
a (l/5;7r, 3/57r) spin excitation was reported, indicating the 
presence of vacancy order and block antiferromagnetic (AF) 



ground state II23L 12411 . For the SC compound, the INS experi- 
ments have revealed several channels of spin excitations |125!- 
I27I1 around the {\l5n, 3/5n), (n, n/2), and (tt, 0), respectively. 
Among them the (1/5:^, 3/5:^) excitation accounts for the va- 
cancy ordered phase. The presence of (7r,0) excitation seems 
controversial: It was indicated in Refs. i25[l26ll that no spin 
excitation at this wave vector was observed, but, on the con- 
trary, the (n, 0) spin excitation was very recently reported in 
Ref. i27l[ . It may come from the localized moments in the SC 
state ll27ll . Interestingly, the (n, k/2) spin excitation was found 
to be rather robust in the SC state possibly due to the FS nest- 
ing ll25l - l27ll . Notably, this excitation is enhanced significantly 
at the frequency 0.014 eV, corresponding to the spin resonance 
at this frequency. The resonant spin excitation is a common 
feature for high-Tc superconductors and thus it may be re- 
lated to superconductivity intimately, providing an insight for 
exploring the mechanism and pairing symmetry of supercon- 
ductivity. 

To study spin excitations and other physical properties the- 
oretically, it is highly-demanded and desirable to establish a 
minimum model that is able to capture the essential features 
of this compound. Firstly, the vacancy order can be neglected 
when studying the physics of superconductivity. It has been 
reported by several groups that the phase separation occurs 
in the SC state, with a vacancy ordered state only existing in 
the insulating region [I28I432I1 . Based on the ARPES experi- 
ments |i4n8i| there is no band folding effect corresponding to 
the vacancy order wave vector, also seems to propose that the 
vacancy order does not exist in the SC state. As a result, the 
vacancy order should not be relevant to superconductivity and 
need not to be concerned for a minimum model. 



From the band calculations for this compound 1191- 11 ID . it 
seems that all of the five orbitals hybridize strongly, and thus 
it was proposed that all five c/-orbitals should be considered 
to construct a model Note that the five-orbital 

model is one of the most frequently-used model to describe 
the iron-pnictide compound ll33ll . However, this model in- 
cludes some non-essential parts that may be redundant for 
superconductivity (especially those bands far away from the 
FS), and even makes it quite difficult to accurately analyze 
and understand certain essential physics behind because too 
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many unknown parameters are involved. For example, if the 
electron correlation needs to be included, both inter-orbital 
and intra-orbital interactions should be taken into account; 
while the corresponding interaction strengths are difficult to 
be determined (or adjusted) in order to figure out the relevant 
physics. Motivated by this, a minimum two-band model was 
proposed for iron pnictides (3^j3?], and many physical prop- 
erties have been understood based on this type of two-band 
model, especially for those related closely to the low-energy 
excitations. When studying the A^Fe2_vSe2 materials, a simi- 
lar two-orbital model was put forward, while the obtained FS 
size is much larger than that obtained from the band calcula- 
tion 1 36]. Intriguingly, as discussed below, we can refine the 
minimum model of this compound to have one band to cap- 
ture the essential physics with better results. As seen from 
the band calculations ll9l-[Tlll. there exist two electron-like FS 
sheets around the M point. Around the F point, a small FS 
pocket may exist, but it is k- dependent and electron-like. The 
FSs from ARPES experiments are consistent with the LDA 
calculations while the electron pockets around the F point 
have a very low spectral weight. So it is reasonable to be- 
lieve that only the FS pockets around M point are essential 
for constructing a minimum model, while all other bands that 
do not cross the Fermi energy may be neglected. Taking into 
account that one unit cell consists of two iron ions plus the 
band folding effect, only one kind of FS pocket around the 
X = (tt, 0) or its symmetric points is actually relevant in the 
unfolded BZ. 



In this paper, motivated by the above considerations, we 
propose a single-band tight-banding model as a minimum one 
for the SC A;cFe2-ySe2. Our model could have a band disper- 
sion crossing the FS, being qualitatively consistent with that of 
the five-orbital model |22, 37] fitted from the LDA band cal- 
culations and the ARPES experiments for the A iFe2-vSe2 ma- 
terials [4- 8]. Thus the present model may serve as an effective 
one for describing the low energy physics of the AfFe2-vSe2. 
The normal state spin susceptibility is calculated and then an- 
alyzed based on the FS topology. The SC gap is also calcu- 
lated self-consistently, with the robust c/-wave SC pairing be- 
ing revealed. The cZ-wave symmetry can well be understood 
by the spin fluctuation picture and the fermiology theory. The 
SC state spin excitations are investigated further with the d- 
wave symmetry. The resonant spin excitation at (tt, tt/I) is ob- 
served, which is consistent with recent INS experiments and 
thus provides a certain kind of support for the cZ-wave symme- 
try in AvFe2-vSe2 SC materials. 



The paper is organized as follows. In Sec. II, we introduce 
the model and work out the formalism. In Sec. Ill, we per- 
form numerical calculations and interpret the obtained results. 
The band structure, pairing symmetry and spin excitations are 
elaborated, respectively. Finally, we give a brief summary in 
Sec. IV. 



II. MODEL AND FORMALISM 

We start from at- J type model including the tight-banding 
term and the local spin interaction, which reads 



H : 



ko- 



H, 



(1) 



where et is phenomenologically taken as the single band 
tight-banding approximation: Sk = -2f(cosA:v + cos A;,,) - 
At' cos k f cos ky - ji. Hj is the local spin interaction. From the 
first principle calculation in Ref. [11] , we assume the spin in- 
teraction to include the nearest-neighbor and the next-nearest- 
neighbor interactions ifTlll : 



Hj^J,Y^SrSj + J2Y^SrSj. 



(2) 



<'7> 



('■;>' 



Here {ij) and {ij)' represent the summation over the nearest 
and next-nearest neighbors, respectively. In the present work, 
we set 7i : 72 = 1 : 0.9. 

The bare spin susceptibility determined from the tight- 
banding part can be calculated as. 



^o(q, oj) 



-y 

N 4^ 



/(gk+q) - /(gk) 

^^-(gk+q-gk) + ir' 



(3) 



where f(x) is the Fermi distribution function. 

The correction of the spin fluctuation induced by the spin- 
spin interaction Hj may be included in the random phase ap- 
proximation (RPA), 



1 + J^o{q,Loy 



(4) 



where Jq = Ji (cos kj, + cos ky) + 2J2 cos kjc cos ky is the Fourier 
factor of the spin coupling term Hj. 

Considering that the SC pairing is meditated by the spin 
fluctuation, we may write the linearized eliashberg's equation. 



^A(k) 



y v(k-k') 



, tanh(j6ek,/2) 



k 



2gk' 



A(k'), 



(5) 



where /? = 1 /T. We consider the spin-fluctuation as the effec- 
tive pairing potential V(q) - g^xi^.^ 0)- Since we here address 
the pairing symmetry, we focus only on the zero-energy spin 
susceptibility, which should produce qualitatively correct re- 
sults for the pairing symmetry as usual |38l'39t]. 

With the SC order parameter, the bare spin susceptibility in 
the SC state can be expressed as, 

1 Y^,lr, gkgk+q + AkAk+q /(^k+q) - /(^k) 

^o(q,w) = - + T-^ ] ^ . j7 . -T 

N ^ 2 £k£k+q W - £k+q +E\, + lT 

1 gkgk+q + AkAk+q 1 - /(^k-q) " /(^k) 

+ -[\ -] 

4 ^k^k+q W + Ek+q +E],+lT 

_ gkgk+q + AkAk+q ^ /(^k+q) + /(^k) " 1 ^ 

4 ^k^k+q - £^k+q " £^k + 

where Ek is the quasiparticle energy in the SC state with 
/ik - -^g^ + A^. The renormalized spin susceptibility in the 
framework of RPA can be obtained by Eq.(4). 
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FIG. 1 : (a) The band structure along different momentum cut in the 
extended BZ. (b) A replot of the band structure in the folded BZ with 
considering one unit cell including two Fe ions. 



III. RESULTS AND DISCUSSIONS 

A. Band structure 

Based on the five-band model, there are three bands below 
the Fermi energy and one band above it ll22ll . It is reasonable 
to concentrate on the band that crosses the FS, which plays 
an essential role in superconductivity. Considering that the 
electron filling per Fe site is 6 + ^, the three filled bands con- 
tribute six electrons per site and thus the electron filling for 
the band that crosses the Fermi energy is 6. These features 
are crucially captured by the the present single-band model. 
We depict in Fig. 1(a) the band dispersion with the hopping 
constants t and t' being set as 0.02 eV and -0.045 eV, re- 
spectively. The chemical potential is controlled by the doping 
density 6 with 6 - 0.2. Considering one unit cell including 
two Fe irons, we plot the band dispersion in the folded bril- 
louin zone in Fig. 1(b). The low energy band dispersion is 
qualitatively consistent with the five-orbital model and LDA 
calculation li22l IbtII . While here the band width is much 
smaller than those obtained from five-orbital model. Actu- 
ally, the five-orbital is usually fitted from the LDA calculation, 
which should generate qualitatively correct FS while usually 
the band width is larger due to neglecting the electron corre- 
lations. In the present work, a smaller renormalized hoping 
constant is used and the obtained band width and doping den- 
sity are consistent with the ARPES experimental results 



B. Pairing symmetry 

The pairing symmetry is one of the most important issues, 
which can provide us the information of the pairing mecha- 
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FIG. 2; (a) The bare and renormalized zero energy spin susceptibility 
with different 7i along different momentum cut. (b) The intensity 
plot of the renormalized spin susceptibility x(<l, = 0) with = 
0.25 ev. 



nism. While this is still under debate for AfFe2-, Se2. Experi- 
mentally it has been revealed that the order parameter is node- 
less and nearly isotropic along the FS pockets 0-01 • Since no 
phase-sensitive measurement has so far been reported, both s- 
wave and nodeless cZ-wave symmetries may be able to repro- 
duce the experimental results. Theoretically, some predicted 
that the order parameter should have d-wave pairing fl3-[T^ 
and some proposed the 5- wave symmetry |18, 19]. There are 
also some theories suggested that the pairing symmetry is ei- 
ther li-wave or .s-wave, or even a mixture of them i20l - t22ll . 

We now study the pairing symmetry based on the spin dy- 
namics scenario. In the RPA framework, the spin suscepti- 
bility are mainly determined from the two factors: one is the 
bare spin susceptibility xo evaluated from the tight-banding 
term and the other is the RPA factor 1 + J^i) closely related 
to the local spin interaction. The behavior of the bare spin 
susceptibility is usually determined by the FS topology, with 
the maximum spin excitation occurring at the nesting wave 
vector For the renormalized one, the pattern of Jq, represent- 
ing the spin interaction strength in the Fourier space, plays an 
important role. While at present the parameters 7, in Jq are 
difficult to be determined theoretically. They usually deviate 
from the actually value of interchange integral constants. In 
cuprates it is estimated roughly from the AF critical point 
or the spin resonance frequency 1.41,1 . In the present work, we 
consider several values of Ji to explore the renormalized spin 
susceptibility. 

The bare and renormaUzed spin susceptibility along the 
two-dimensional cut are plotted in Fig. 2(a). The intensity 
plot of the renormalized spin susceptibility in the whole BZ 
is plotted in Fig. 2(b). As is seen, the bare spin susceptibility 
shows two broad features around the momentums (0, 0) and 
(tt, n). Here the maximum spin susceptibilities at four wave 
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FIG. 3: (a) The gap function from the self-consistent calculation with 
the coupling strength g = 0. 1 eV. (b) The gap along the FS sheets 
with the dashed and solid lines are around the X and Y points, re- 
spectively, with denoted in panel (a). 



vectors Qi_4 are revealed. In the RPA approach, the spin ex- 
citations around (0, 0) are suppressed and those around (;r, n) 
enhanced, respectively. As Ji increases to 0.25 eV, the maxi- 
mum spin excitation appears at the momentum Q3 = (n, n/T). 
The (tt, 7r/2)-spin excitation can be seen more clearly from 
Fig. 2(b). As displayed, four peaks appear at the wave vector 
{7T,n/2) and its symmetric points. This result is qualitatively 
consistent with the recent INS experiments ll25il26ll . A further 
comparison between our theoretical results with INS experi- 
ments will be presented in Sec.III(C). 

We here look into the pairing symmetry from Eq.(5), 
namely, the temperature T is decreased until the maximum 
eigenvalue /I = 1 is obtained, with the SC gap being the eigen- 
vector for the maximum eigenvalue. The intensity plot of the 
SC gap is displayed in Fig. 3(a). Obviously, the pairing sym- 
metry is of t/j.2_^2-wave. We have checked numerically that 
this result is rather robust to the reasonable parameter change. 
The SC gap along the FS for the d-wave pairing symmetry 
A - A()/2(coskx-cosky) is plotted in Fig. 3(b). As is seen, the 
SC gap along one FS sheet is nearly isotropic with the differ- 
ence less than 10%, being consistent with the ARPES exper- 
iments ||3-13l- In addition, it was reported in ARPES experi- 
ments that the gap is small or vanishing for the smaller pockets 
around the F point |4, 5], being also qualitatively consistent 
with the (i^2_,,2-wave. However, it was indicated in Refs. IHSl 
that there exists a larger FS pockets around F with the gap 
being nearly isotropic, which seems to contradict with the d- 
wave symmetry. While we think that the larger FS around F 
may come from the band folding effect: dominant features of 
the pocket and the gap may be folded from those around M 
point. If this is the case, our results for the gap symmetry 
agree qualitatively with the ARPES experiments. 

The above results of the spin fluctuation can be addressed 
soundly based on the fermiology picture. The FS is plotted in 
Fig. 4. The FS sheets around X and its symmetric points are 
clearly seen. We denoted the nesting wave vectors (Q1-4) in 
Fig.4 |42]. As seen, all of the maximum bare spin excitations 
shown in Fig. 2a are related to the FS nesting. Based on the 
fermiology and Eq.(3), it is rather clear that the spin suscep- 
tibility reaches its largest value at Qj because fik+Q, - Sk is 
vanishingly small. 

We now elucidate the origin of the d-wave symmetry. At 
the RPA approach, the spin excitations at Qi,2 from intra- 
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FIG. 4: (a) The normal state FS with the arrows denoting the different 
nesting vectors. 



pocket FS nestings are suppressed and those at Q3_4 from the 
inter-pocket ones enhanced. Thus the inter-pocket scattering 
should play major role to achieve superconductivity. The pair- 
ing potential V contributed by the spin fluctuation is largest at 
the wave vector Q - Qxi- The factor tanh(y6ek'/2)/2ek' in 
Eq.(5) is positive for any Sk' and largest at £k' = 0, which 
means that the pairing near the FS is important. The gap func- 
tion at or near the FS should satisfy the condition Ak = -Ak+Q 
according to Eq.(5). As shown in Fig.4, if k belongs to one 
sheet of FS, then k + Q should be near the other neighbor- 
ing sheet of FS. For the t/(.2_y2-symmetry, as displayed in 
Fig.4(b), the SC gaps have the same magnitudes and differ- 
ent signs along the two neighboring FS sheets. The condition 
Ak - -Ak+Q is satisfied approximately. In this sense, we give 
an intuitive understanding of the pairing symmetry in this ma- 
terial. 



C. Spin excitations in tlie SC state 

Now let us study the spin excitations in the SC state. Here 
the d-wave SC order parameter Ak - A()/2(cosk^ - cos ky) 
with the maximum gap Aq = 0.01 eV being considered. The 
imaginary parts of the spin susceptibility Im;^'(q, w) as a func- 
tion of the wave vector q are plotted in Fig. 5. One can see 
clearly the maximum spin excitation near the wave vector 
Q = (tt, 7r/2) [with the exact peak center at (tt, 0.5 IStt)] at 
o) - 0.014 eV. As the frequency increases to co - 0.018 meV, 
the intensity of the spin excitation is enhanced. The maxi- 
mum excitation persists at the wave vector Q [The exact peak 
center is at (tt, 0.5357r)]. The spin excitations near the wave 
vector (tt, 7r/2) is well consistent with the experimental re- 
sults revealed by the INS experiments 1I25I l26ll . Moreover, 
the slight shift of the peak center as the frequency increases 
and the weak dispersion behavior are also consistent with the 
experimental results [26]. 

In the following, we study the frequency dependence of the 
spin excitations at the wave vector Q = (tt, 7r/2). Here three 
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FIG. 6: (a) The imaginary parts of the spin susceptibility versus the 
energy for the wave vector Q = (n,n/2) in the normal state and 
SC states with different symmetries, respectively. (b)The bare spin 
susceptibilities versus the energy. The black line denotes the real 
part and the red line imaginary part, respectively. The green line is 
the real part of the RPA factor. 



FIG. 5: The intensity plots of the imaginary part of the spin suscep- 
tibility Im;^' as a function of the momentum with the energies 0.014 
eV and 0, 018 eV, respectively. 



kinds of pairing symmetry, namely, the t/j.:_,,2-wave pairing 
symmetry with Ak = Ao/2(cosA;j. - cos A:,,), 5 -pairing sym- 
metry with Ak = A() cos kx cos ky, and isotropic i-wave pairing 
symmetries with Ak = Ao, are considered. As is known, all 
of the three pairings would generate isotropic SC magnitude 
along the FS. The imaginary parts of spin susceptibility Im;^' 
as a function of frequency in the normal state and SC state 
are presented in Fig. 6(a). We observe clearly that for the 
(i-wave pairing symmetry, the spin susceptibility at the fre- 
quency 0.018 eV is enhanced significantly, indicating the spin 
resonance for this energy. For the and isotropic i-wave 
pairing symmetry, the low energy spin susceptibilities are al- 
ways less than those in the normal state. Thus there is no reso- 
nant spin excitation for these two symmetries. This significant 
difference between the li-wave symmetry and .s-wave symme- 
try is interesting and this feature can be used to determine the 
pairing symmetry through comparing with the experiments. 

The above difference of the spin excitations can be under- 
stood through analyzing the coherence factor C in Eq.(6) with 
C - 1 - "''^''^-''t.^''^''^'' . The spin excitation is enhanced as 
AkAk+Q < (with k and k -i- Q being the momentum close 
to the FS). This condition can only be satisfied for the c/-wave 
pairing symmetry, which provides a natural explanation for 
the resonant spin excitation. The origin of the spin resonance 
can be clarified further in the framework of RPA. The bare 
spin susceptibility and the PRA factor with the d-wave pair- 



ing are plotted in Fig. 6(b). As seen, due to the presence of 
the SC gap, the imaginary part of the bare spin susceptibility 
approaches to zero at low energies, ascribed to the spin gap. 
At the edge of the spin gap (near 2A()), it has a steplike rise. 
In the mean time, the real part of the bare spin susceptibility 
ReXo develops a sharp structure and reaches the maximum at 
this frequency. Thus the RPA factor reaches the minimum (or 
sometimes equals to zero coiTesponding to a strong resonant 
state at the frequency oj < 2A()) at this frequency, which plays 
the major role for the appearance of the spin resonance state. 

Recently the spin resonance at the wave vector Q = (n, 7t/2) 
and the frequency 0.014 meV is observed by the INS experi- 
ments 1 25- 27]. Our numerical results for the c/-wave pairing 
symmetry are qualitatively consistent with the experimental 
results, while those for the .s-wave pairing symmetry contra- 
dict obviously with the experiments. In this sense, we have 
provided likely an indication that supports the li-wave pairing 
symmetry in this family of SC materials. 



IV. SUMMARY 

In summary, we have established an effective single-band 
model, which captures likely the essential low-energy physics 
in the AiFe2-vSe2 material. Several puzzled properties have 
been explained satisfactorily based on this minimum model. 
An intriguing spin excitation with the wave vector {n,jil2) has 
been revealed. In addition, from the fermiology analysis, we 
have developed a coherent picture for the spin excitations and 
the unconventional li.iJ-v- pairing symmetry. The spin exci- 
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tations in the SC states are studied and the spin resonance at 
{tt, njl) is revealed for the d^2_^2 symmetry. While there is no 
resonant spin excitations for the s-wave symmetry. Thus we 
have provided likely an indication that supports the t/-wave 
symmetry in Aj:Fe2-ySe2 materials. All of our results are in 
qualitative agreement with the ARPES and INS experiments. 
In addition, the simple effective model presented here is quite 
useful and promising for exploring rich but non-trivial physics 
in iron-selenide systems. 
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